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Context. The cosmological 21 cm line promises to be a formidable tool for cosmology, allowing the investigation of the end of the so-called 
dark ages, when the first galaxies formed. 

Aims. Astrophysical foregrounds are expected to be about three orders of magnitude greater than the cosmological signal and therefore represent 
I a serious contamination of the cosmological 21 cm line. Detailed knowledge of both their intensity and polarization structure on the relevant 
, angular scale of 1-30 arcmin will be essential for extracting the cosmological signal from the data. 

Methods. We present the first results from a series of observations conducted with the Westerbork telescope in the 140-160 MHz range with a 
, 2 arcmin resolution aimed at characterizing the properties of the foregrounds for epoch of reionization experiments. The polarization data were 
' analysed through the rotation measure synthesis technique. We computed total intensity and polarization angular power spectra. 
Results. For the first time we have detected fluctuations in the Galactic diffuse emission on scales greater than 13 arcmin at 150 MHz, in the 
low Galactic latitude area known as Fan region, centred at o- = 3*' 10'", 6 = 65°30'. Those fluctuations have an rms of 14 K. The total intensity 
power spectrum shows a power-law behaviour down to ~ 900 with slope yS' = -2. 2 ±0.3. The detection of diffuse emission at smaller angular 
scales is limited by residual point sources. We measured an rms confusion noise of ~3 mjy beam"'. 

Diffuse polarized emission was also detected for the first time at this frequency. The polarized signal shows complex structure both spatially 
and along the line of sight. The polarization power spectrum is not affected by residual point sources and is only limited by the thermal noise. 
It shows a power-law behaviour down to C ~ 2700 with slope = -1.65 ± 0.15. The rms of polarization fluctuations is 7.2 K on 4 arcmin 
. scales. 

' Conclusions. The measured total intensity fluctuations are used to estimate the foreground contamination on the cosmological signal. By 



extrapolating the spectrum of total intensity emission, we find a contamination of ST = -^({f + l)C^/2;r ~ 5.7 K on 5 arcmin scales and a 
corresponding rms value of ~18.3 K at the same angular scale. 

The level of the polarization power spectrum is 6T ~ 3.3 K on 5 arcmin scales. However, the Fan region cannot be taken as representative of 
other sky regions, given its exceptionally bright polarized signal, but is likely to represent an upper limit on the sky brightness at moderate and 
high Galactic latitude. 

Key words. Polarization - Cosmology: diffuse radiation - Cosmology: observations - Radio continuum: general - ISM: general - ISM: 
magnetic fields 



1 . Introduction the cosmic microwave background polarization suggest that a 
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ABSTRACT 




Although interesting in its own right, the study of the dif- 
fuse Galactic radio emission at metre wavelengths has recently 
taken on more significance because this emission is the dom- 
inant contaminating foreground for the study of the redshifted 
21 cm line from the epoch of reionization (EoR, Shaver et al. 
1999; Furlanetto et al. 2006). The recent WMAP results on 



significant part of the reionization of the Universe must have 
occurred at z < 13, while the Gunn-Peterson effect in high- 
redshift quasars tells that the bulk of the reionization must have 
occuiTed earlier than z~6. The best frequencies for search- 
ing for the cosmological 21 cm line therefore lie in the 100- 
200 MHz range. 
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The 21 cm line is expected to have an rms brightness tem- 
perature of a few millikelvin on angular scales of a few arcmin 
(Zaldarriaga et al. 2004). However, both the angular scale and 
the intensity depend on the details of the reionization process. 
Simulations indicate that the cosmological signal is expected to 
vary with redshift (Ciardi & Madau 2003; Mellema et al. 2006; 
Thomas et al. 2008). 

The Giant Metrewave Radio Telescope (GMRT) and sev- 
eral dedicated experiments like LOFAfQ, MWaH and 21CMA0 
will search for these signals. 

The detection of cosmological signals with millikelvin 
brightness temperatures is plagued by a wide range of chal- 
lenges. Foremost among these is the intense foreground emis- 
sion from our Galaxy (Shaver et al. 1999). In fact about 75% 
of the sky brightness is due to diffuse emission from our 
Galaxy, most of which is caused by the synchrotron process. 
Our present knowledge of the Galactic diffuse radio emission 
at metre wavelength is not adequate to estimate the contami- 
nating effects on the EoR signal. Several all-sky maps (Haslam 
et al. 1982; Reich & Reich 1986; Reich et al. 2001; Page et al. 
2007) at different frequencies and angular resolutions are avail- 
able, but the radio emission from our Galaxy at the relevant 
frequencies (100-200 MHz) is poorly known. The 150 MHz 
single-dish all-sky map by Landecker & Wielebinski (1970) 
has only 2° resolution, while estimates of the spatial fluctua- 
tions on arcmin scales are required in EoR experiments. 

The first indications that the Galactic foreground lacks sig- 
nificant fine-scale structure at low frequencies were presented 
by Wieringa et al. (1993). Subsequent dedicated arcmin reso- 
lution observations with the Westerbork telescope (WSRT) at 
frequencies from 315 to 380 MHz, have indeed shown that the 
diffuse synchrotron emission from the Galaxy is generally very 
smooth and therefore mostly resolved out by the interferom- 
eter. This also became evident from the Westerbork Northern 
Sky Survey (WENSS) at 325 MHz (Rengelink et al. 1997, see 
also Wo WO). 

Extrapolating these limits to lower frequencies and smaller 
angular scales is uncertain because both the spectral and the 
spatial properties of the diffuse emission can vary across the 
sky. The mean spectral index of the synchrotron emission at 
high Galactic latitude has been quite well constrained to be yS = 
2.5 +0. 1 in the 100-200 MHz range (Rogers & Bowman 2008). 

The spatial properties are less well known because no 
power spectra of the diffuse emission exist on arcmin scales. 
La Porta et al. (2008) analysed power spectra of the diffuse 
Galactic synchrotron radiation down to scales of a few degrees 
and found significant variations across the sky even after sub- 
tracting the brightest sources. 

Although the redshifted 2 1 cm line signal is itself unpolar- 
ized, we still have to worry about intrinsically polarized sig- 
nals from the sky. The instrumental polarization of low fre- 
quency radio telescopes and the upcoming arrays is significant 
(10-20%). Imperfectly calibrated polarized signals from the 

' http://www.lofar.org 

^ http://haystack.mit.edu/ast/arrays/mwa 
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* http:/www.astron.nl/~wow 



sky will then leak into the total intensity images. Even small 
variations in the rotation measure of the polarized signals wiU 
then lead to frequency dependent residuals, corrupting the fee- 
ble cosmological 21 cm signal. It is therefore prudent to search 
for EoR target fields with low foreground polarization. 

Until recently, most studies of Galactic polarization were 
done at relatively high frequencies. Brouw & Spoelstra (1976) 
collected data with the Dwingeloo telescope, which cover large 
regions of the sky in the 408-141 1 MHz range with a resolution 
between 0.5° and 2°. A full-sky absolutely-calibrated polar- 
ization map at 1.4 GHz with 0.5° resolution was recently pub- 
lished by Wolleben et al. (2006). Selected regions of low total 
intensity synchrotron emission have been observed between 1 .4 
and 2.3 GHz (Bernardi et al. 2003; Carretti et al. 2005a; Carretti 
et al. 2005b; Carretti et al. 2006; Bernardi et al. 2006) charac- 
terizing the spatial properties of the diffuse polarized emission 
down to arcmin scales with the purpose of studying the con- 
tamination of the cosmic microwave background polarization. 

At lower frequencies, the sky coverage has been scarce. 
Following the pioneering WSRT work by Wieringa et al. 
(1993) at 325 MHz, several dedicated 350 MHz observations in 
selected regions have been conducted (Haverkorn et al. 2003a; 
Haverkorn et al. 2003b - hereafter HKB) together with observa- 
tions of a limited number of WENSS survey fields (Schnitzeler 
2008). These studies indicated that foreground polarization lev- 
els, even at higher Galactic latitudes, could easily go up to 5 K 
or more on 5-10 arcmin scales (de Bruyn et al. 2006). If these 
regions were Faraday thin, the polarized brightness tempera- 
tures at 150 MHz could go up to many tens of kelvin. This 
would place extraordinary requirements on the level of accu- 
racy of the instrumental polarization calibration. 

Several authors have simulated the properties of the fore- 
grounds in order to test strategies and algorithms to subtract 
or filter them from the data (Di Matteo et al. 2004; Morales et 
al. 2006; Gleser et al. 2007; Jelic et al. 2008, Liu et al. 2009a, 
Harker et al. 2009, Liu et al. 2009b). It is clear, however, that 
these simulations require guidance from observational data on 
both total and polarized foreground levels. 

We have therefore initiated a comprehensive programme to 
directly measure the properties of the Galactic foreground in 
the 100-200 MHz range, the relevant frequencies for the EoR. 
We used the Low Frequency Front Ends (henceforth LFFE) 
on the WSRT to carry out these observations. We chose three 
different fields. One field is centred around the so-called Fan 
region, at Galactic coordinates I - 137° and b = +8° in 
the 2nd Galactic quadrant. It was selected because it remains 
highly polarized down to the 350 MHz frequencies (Brouw & 
Spoelstra 1976; Haverkorn et al. 2003a). The other two fields 
were chosen because they represent possible targets for EoR 
observations. The second field is centred on the very bright ra- 
dio quasar 3C196 in one of the coldest regions of the Galactic 
halo. The third field is located close to the North Celestial 
Pole, which would allow us to collect night time observations 
throughout the year at the geographic latitude of the LOFAR 
an-ay (h-53°). 

In this paper we present results from WSRT observations 
of the Fan region at 150 MHz aimed at studying the proper- 
ties of both the total intensity and the polarized diffuse emis- 
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Table 1. Summary of the observational setup. 



Coordinates of the field 

centre (J2000.0) 

Number of spectral bands 

Central frequency of each band (MHz) 

Width of each band (MHz) 
Frequency resolution (kHz) 
Time resolution (sec) 
Angular resolution 
Conversion factor 



a = 340", 5 = 65°30' 
8 

139.3, 141.5, 143.7, 145.9, 
148.1, 150.3, 152.5, 154.7 
2.5 

4.9 (9.8 after tapering) 
10 

2' X 2'cosec(^) ^ 2' x 2.2' 
1 mjy beam ' = 3.98 K 



sion as contaminants for the EoR signal. In subsequent papers 
we will present results from the analysis of the other fields, 
together with the challenges and the calibration strategies for 
foreground subtraction applied to the WSRT data. 

This paper is organized as follows: in Section|2]we present 
the observations together with the data reduction and the basic 
results, in Section [3] we describe the results of the total inten- 
sity emission, in Section |4] we analyse the polarized emission, 
in Section |5] we present the power spectrum analysis and in 
Section|6]we conclude. 

2. Observations and data reduction 

Observations of the Fan region took place between the end of 
November and the beginning of December 2007 for a total of 
6 X 12 hours. All observations of the target field took place 
between sunset and dawn, avoiding coiTuption of the shortest 
baselines by the strong and often time-variable emission from 
the Sun. Short calibration observations of the highly polarized 
pulsar J0218+4232 (Navarro et al. 1995) and the strong unpo- 
larized source 3C196 were made immediately before and after 
the target observation, which proceeded uninterrupted. 

The WSRT telescope consists of 14 dishes of 25 m diame- 
ter, ten of which (labelled to 9) are on fixed locations 144 m 
apart. The other four (labelled A to D) are movable on a rail 
track. The redundant baselines are normally not included in the 
imaging process, because they cause very strong grating lobes. 
The normal choice is to use all the correlations between the 
ten fixed antennas and the four movable ones. Using six con- 
figurations of the movable telescopes, the four movable dishes 
were incrementally moved by 12 m obtaining a uniform uv 
coverage from 36 m up to the maximum spacing of 2760 m. 
Unfortunately antenna 5 was missing in all the observations, 
causing two broad gaps in the uv plane centred at ~ 300 |u| and 
~ 900 |u|, where |u| = Vm^ + is the projected distance in the 
uv plane measured in wavelengths. 

Table [T] summarizes the main characteristics of the obser- 
vations. The eight spectral windows were chosen to provide a 
contiguous frequency coverage, allowing for a small overlap 
between the bands. A large part of the fourth spectral band, 
centred at 145.9 MHz, turned out to be too severely contami- 
nated by radio interference and the whole band was discarded 
from the analysis. 

After Manning tapering the 512 channels, and discarding 
the edges of each band, we were left with a total of about 



240 frequency channels per band. Each channel has a width 
.of 9.8 kHz. 

The data were reduced using the AIPS-H-tH package. We 
integrated the standard AIPS-I-+ distribution with flagging rou- 
tines explicitly developed for dealing with radio frequency in- 
terference (RFI) and with routines which determine the polar- 
ization calibration. Dipole gain and leakage corrections were 
determined using the unpolarized calibrator 3C196. Due to an 
unknown phase offset between the horizontal and the vertical 
dipoles, signal can leak from Stokes U into Stokes V (Sault, 
Hamaker & Bregman 1996). In order to correct for this phase 
'diff'erence, the polarized pulsar PSRJ02 18 -1-4232 was observed 
for about 15 min just prior to the start of the 12 hour syn- 
thesis. Since the pulsar has a known rotation measure RM - 
-61 rad mT^ (NavaiTO et al. 1995), we corrected the phase dif- 
ference by rotating the polarization vector in the plane defined 
by the Stokes U -V parameters in order to have zero Stokes V 
flux, where the direction of the rotation has to provide a nega- 
tive RM for the pulsar 

The flux scale at low frequencies is not very accurately de- 
termined and a more definitive flux scale determination is be- 
ing developed for LOFAR. In our case we used the radio source 
3C196 as the primary WSRT flux calibrator at low frequencies. 
3C196 is a very bright steep spectrum radio galaxy whose spec- 
trum appears smooth down to about 10 MHz (Laing & Peacock 
1980). We have adopted a flux density at 150 MHz of 76.8 Jy 
for 3C196 and a power-law spectral index of a - -0.64 in the 
frequency range from 1 15-175 MHz - although there may be a 
slight steepening at the upper end of this band. We believe this 
number to be accurate to ~5%. 

The system noise is constituted by the receiver noise and 
the brightness temperature of the sky itself. At low frequency 
our Galaxy contributes significantly to the system noise. 3C196 
is located in the coldest spot of the Galactic halo according to 
the 408 MHz all-sky map (Haslam et al. 1982). The System 
Equivalent Flux Density (SEED) of the telescopes at 150 MHz 
is about 7000 Jy in the 3C196 area but it rises to over 10,000 Jy 
in the Galactic plane where the Fan region is located. Because 
the WSRT receivers operate with an automatic gain control 
system before the analog-to-digital converter, it continuously 
measures the total power to allow corrections for the variable 
input levels. 

Since the total power detectors (which integrate the power 
over the whole 2.5 MHz band) are corrupted by RFI for most 
of the time, we cannot automatically correct the correlation co- 
efficients for the variations in the system noise. The LFFE band 
is full of mostly impulsive and narrow band RFI coming from 
airplanes, satellites and mobile users as well as electronic hard- 
ware within the building which is located halfway along the ar- 
ray. At the high spectral resolution of 9.8 kHz provided by the 
backend most of this RFI can be excised. However, the total 
power data must be manually inspected for suitable stretches 
of power level measurements. These data form the basis for a 
manual correction of the flux scale. 

We found that the total power ratio between the FAN region 
and the 3C196 field is 1 .48 ± 0.03 over the 138-157 MHz band. 



http;//aips2. nrao.edu/docs/aips -H+.html 
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Following the transfer of the bandpass correction determined 
for 3C196 we have therefore applied an additional correction 
of a factor 1.5 to the visibility data. 

After this correction was applied, the data were self- 
caUbrated in order to correct for time variations in the dipole 
gains. Self-calibration is especially useful at metre wavelength 
because radio interference causes temporal variations in the 
gains and ionospheric turbulence introduces temporal varia- 
tions in the visibility phases. 

The signal to noise ratio (SNR) in our data is quite poor, 
given that the brightest source is ~2.8 Jy and the noise per visi- 
bility, per channel, per time slot and per polarization is ~30 Jy. 
Moreover we want to avoid averaging the visibilities in time to 
improve the SNR because a time average decorrelates the visi- 
bilitity signal if the ionosphere is turbulent on short time scales. 
We therefore adopted the following procedure to calibrate the 
data. 

After the bandpass correction determined from 3C196 was 
applied, we generated frequency cube images for every night 
and every spectral band. We normalized the dirty images to the 
average peak value derived from every cube. In this way we 
compensated for frequency dependent amplitude errors. 

Cyg A and Cas A were then subtracted from the data be- 
cause at these frequencies they are bright enough to generate 
significant side lobes within the field of view - at a level of 
~150 mJy for Cas A - despite being far outside the target field. 
Afterwards we created a sky model per night and per spectral 
window which contained the brightest 70-80 sources, down to 
a flux limit of ~300 mJy. For every night and for every spec- 
tral window we then averaged all the channels to improve the 
SNR, and computed phase solutions every two minutes to cor- 
rect for ionospheric turbulence. After this phase correction was 
applied, the data in the spectral bands centred at 143.7 and 
148.1 MHz still showed amplitude variations as a function of 
time. We corrected for amplitude variations in the complex 
gains by computing solutions every five minutes. 

The resulting Stokes / map obtained by averaging aU the 
seven spectral windows is shown in Figure [T] Since the Half 
Power Beam Width (HPBW) of the WSRT dish is -6° at 
150 MHz, Figure [U maps the full primary beam. Unless stated 
otherwise, none of the the images shown in this paper are cor- 
rected for the primary beam attenuation. The noise is therefore 
constant across the images, but the celestial source brightness 
is modulated by the radial dependence of the beam gain. 

The image is mostly dominated by emission from unre- 
solved point sources. Faint diffuse emission is visible in the 
South- West corner of the image where the Galactic Plane is lo- 
cated. The brightest emission regions, which are heavily atten- 
uated by the primary beam, are due to HII regions and diffuse 
structures within the Perseus arm located near Galactic coordi- 
nates ^7 ~ 1° and 133° < I < 138°. 

3. Diffuse Galactic foreground 

Since our main interest is to characterize the fluctuations 
in the diffuse Galactic foreground, we removed the point 
sources from the map. Several thousand sources were identi- 
fied through a CLEAN deconvolution in a 12° x 12° image. 



The CLEAN model, containing all the sources brighter than 
15 mJy, was then subtracted from the visibility data. Figure |2] 
shows the residual image made from the sum of seven bands. 
The diffuse structure in the Perseus arm now becomes clearly 
visible. The nearby galaxy IC342 is also noticeable as a bright 
structure around a = 3*'38'" and 6 = H-67°30'. 

Figure 12] still shows error patterns such as rings of negative 
values around the brightest sources, indicating that the calibra- 
tion is not yet perfect on them. Moreover, spiky error patterns 
are also present around bright sources, which may be due to 
non-isoplanaticity errors caused by ionospheric phase fluctua- 
tions. A direction-dependent calibration will be applied in the 
future to try to correct for these errors. 

We can measure the noise from the image itself by taking 
the nns value at the edge of the field, where the primary beam 
attenuation has reduced the sky brightness substantially. The 
observed noise at the edge of the image is 0.75 mJy beam"', 
making these observations the deepest available to date at these 
frequencies. 

In the map of Figure |2] diffiise structure on scales of 15- 
30 arcmin is clearly visible. However, it is not clear whether 
diffuse structure exists also on arcmin scales. In order to in- 
vestigate this point, we made an image which only contains 
the high spatial frequencies by selecting the baselines with 
|u| > 150 wavelengths for each spectral band. 

The result is shown in Figure|3]where all the angular scales 
greater than ~12 arcmin are absent: the diffuse structure ap- 
pears almost completely filtered out and the most relevant con- 
tribution left is from point sources. This indicates that the dif- 
fuse structure appearing in the map of Figure |2] is on scales 
greater than 10-12 arcmin and no diffuse power is present on 
smaller scales. 

We estimated the rms due to residual point sources by sub- 
tracting from the high resolution image of Figure [3] all the 
sources down to 15 mJy. We considered the inner 6° x 6° area 
corrected for the WSRT primary beam A(/, y) which is approx- 
imated as: 

A(/,y) = cos^(0.065/r) (1) 

where / is the observing frequency in MHz and y is the angular 
distance from the pointing centre in radians. 

Figure|4]shows the resulting image where the point sources 
have become indistinguishable from the background. This map 
has an nns value o-p^ - 7.2 ± 0.3 mJy beam"'. By taking 
the inner 3° x 3° area where the primary beam correction is 
still small, we determined an rms value of ~4.5 mJy beam"'. 
However, if we carefully select regions which appear to be 
free of instrumental or ionospheric residuals we find an nns 
of ~3 mJy beam"'. 

This number is consistent with the expected confusion 
noise that can be estimated from source counts at higher fre- 
quencies. Deep WSRT observations at 1.4 GHz of the Hubble 
Deep Field measured a confusion noise of 5 //Jy beam"' 
(Garrett et al. 2000). By assuming an average spectral slope of 
a - -0.8 for the sources, the confusion noise cr^^'j expected at 
150 MHz can be predicted from the confusion noise at 1 .4 GHz 
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Fig. 1. Stokes / map obtained by averaging all the data. The colour scale saturates (black) at 75 mJy. The tick marks indicate 
lines of constant Galactic latitude and longitude 



as: 

-^l°f = -L^°(^P~3mJybeam-'. (2) 

We can therefore conclude that the map of Figure|2]is confusion 
limited towards its centre rather than limited by thermal noise. 

Afterwards, we characterized the diffuse emission by ap- 
plying a uv taper to degrade the resolution to ~13 arcmin for 
each of the 7 bands after all the sources down to 15 mJy were 
removed. The resulting image is shown in Figure |5] 



In the 13 arcmin image the contribution from individual 
point sources has completely disappeared whereas the pattern 
of fluctuations dominates the whole image. This pattern does 
not correlate with the distribution of point sources shown in 
Figure [3] indicating that possible unsubtracted sources do not 
contribute significantly to the diffuse emission. 

Considering that the thermal noise is reduced to 
less than 1 K and that the conversion factor is now 
1 Jy beam"' ~ 105.6 K, the hot spots in the image are more 
than ten times above the noise. We conclude that the observed 
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Fig. 2. Stokes / map where the point sources were removed down to a level of 15 mJy level. The conversion factor is 1 mJy beam ' 
= 3.98 K. 



fluctuations represent structure in the Galactic foreground radio 
emission. 

We characterize the fluctuations through their rms value 
computed in the inner 6° x 6° region of the map, after correct- 
ing for the power pattern of the primary beam and avoiding the 
local bright features which are definitely not representative of 
the diffiise foreground emission. In Table |2] we report the rms 
values of the residual maps at full and low resolution. Errors 
on the rms values are computed by splitting the field into two 
halves and taking the average of the differences between these 
two values and the rms of the whole field. 

It is important to note that the rms value on 13 arcmin ac- 
counts for difiiise emission only whereas the value at 2 arcmin 
has a significant contribution from unsubtracted point sources. 

Finally, we analysed the behaviour of the diffuse emis- 
sion as a function of the Galactic latitude. In all Galactic ra- 
dio surveys, the amplitude of the fluctuations decreases at high 
Galactic latitudes. Since our data span several degrees just out- 



Table 2. Rms fluctuations in the Fan field as a function of the 
angular resolution. 

Angular resolution (arcmin) rms value (K) rms value (mJy beam"' ) 



2 
13 



30.5: 
14: 



0.9 
1 



7.7 ± 
128 



0.3 

t9 



side the Galactic plane, we computed the rms as a function of 
Galactic latitude. The result is shown in Figure|6] 

There is a clear drop in the signal as we move away from the 
plane. The rms decreases by a factor of ~2.3 over 6° of Galactic 
latitude. In Figure |7] we plot the same behaviour derived from 
the Landecker & Wielebinski map at 150 MHz which has a 2° 
resolution. We considered a strip 12° wide in Galactic longi- 
tude, centred on the Fan region, and we computed the rms as a 
function of Galactic latitude. The rms drops by a factor of ten 
between low and moderate Galactic latitudes {b > 30°), stay- 
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Fig. 3. Stokes / map made by including only the baselines with |u| > 150 wavelengths. The map has been restored after a CLEAN 
deconvolution. The conversion factor is 1 mJy beam"^ = 3.98 K. 



ing somewhat constant at higher Galactic latitude. We expect 
that the sub-degree fluctuations that we have observed follow a 
similar decrease at moderate and high Galactic latitudes. 

It is therefore relevant to notice that the same region of sky 
in the map by Landecker & Wielebinski shows rms fluctuations 
of ~ 90 K on 2° scales. Since their data come from single dish 
observations, they sample all the angular scales down to the 
beam resolution. The fact that our interferometric data filters 
out most of the large scale power causes the rms to drop by a 
factor of ~6.5 from 2° to 13 arcmin scales. 



4. Polarized emission 

As we mentioned in the introduction, the Fan region is known 

to be polarized from its first observations made with the 
Dwingeloo telescope. Those data already showed a certain uni- 
formity in the polarization angles in this area. HKB observed 



the Fan region at 350 MHz and found significant poliirized 
emission on scales of 5 arcmin. 

We analysed the polarization data at 2 m wavelength 
through the rotation measure synthesis technique (Brentjens 
& de Bruyn 2005). The RM synthesis technique takes advan- 
tage of the Fourier relationship which exists between P{^) and 

/+00 
F(0)e^''*'*'d0 (3) 
oo 

where P{y^) is the complex polarized surface brightness, W{/^) 
is a weighting function and F(0) is the complex polarized 
surface brightness per unit of Faraday depth, A the observing 
wavelength and ^ is the Faraday depth (see Brentjens & de 
Bruyn 2005 for a full description). The RM synthesis performs 
a transformation from the space to the Faraday depth space. 
The output of the RM synthesis analysis is a cube of polarized 
maps at selected values of Faraday depth. The Fourier trans- 
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Fig. 5. Stokes / map where the point sources were removed down to a level of 15 mJy and the uv plane was tapered to achieve an 
angular resolution of 13 arcmin. The conversion factor is 1 Jy beam ' = 105.6 K. 



form of W(A^) gives the RM spread function (RMSF), which 
is the resolution in Faraday depth and depends upon the differ- 
ence in between the two furthermost channels. The RMSF 
is the one-dimensional analogue of the Point Spread Function 
in traditional interferometry. 

Table [3] summarizes the most relevant parameters for the 
polarization analysis. The first three parameters reported in 
Table|3]are defined by the relationships (61), (62) and (63), of 
Brentjens & de Bruyn (2005): 6<p is the resolution in Faraday 
depth, <p^^^ is the maximum scale of Faraday depth to which 



our observations are sensitive and ||0maxll is the maximum RM 
value measurable. 

It is worth noting that our observations are only sensitive to 
Faraday depths which are smaller than the RMSF width, so that 
we lack sensitivity to structures which are extended in Faraday 
depth and we are able to sample only Faraday thin regions, with 
an extension in RM space of ~1 rad m"^. 

Although we corrected for the on-axis instrumental polar- 
ization, the WSRT has a strong off'-axis instrumental polarized 
response which increases with distance from the centre of the 
image and can reach up to 20-30% at the beam half-power 
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Fig. 4. 6° X 6° Stokes / map made by including only the base- 
lines with |u| > 150 wavelengths and where the point sources 
were subtracted down to a level of 15 mJy. The conversion fac- 
tor is 1 mJy beam ' = 3.98 K. The image has been divided by 
the primary beam according to Eq[T] 
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Fig. 6. Rms fluctuations at 13 arcmin as a function of Galactic 
latitude, measured in bins of width 1° of Galactic latitude. 

Table 3. Summary of RM synthesis cube parameters 



6(p 


3 rad m ^ 


Vmax 


0.85 rad m^^ 


ll'/'maxll 


2650 rad m - 


Angular resolution 


4.2 arcmin 


Conversion factor 


1 mJy beam ' = 1 K 



radius. The off-axis polarization is also strongly frequency de- 
pendent and the polarization beam pattern is quite complicated 
at these frequencies. Our sky model for polarization analysis is 
therefore somewhat different from the total intensity model. 




1 1 I I I I 

5 25 45 55 85 
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Fig. 7. Rms fluctuations as a function of the Galactic latitude 
measured from the map by Landecker & Wielebinski (1970). 
The data were grouped into bins of size 15° of Galactic latitude 
and 12° of Galactic longitude. 



We identified all the point sources down to 150 mJy in each 
single channel image and fitted their flux under the assump- 
tion that they were instrumentally polarized, i.e. they did not 
have intrinsic polarization. These sources were then subtracted 
from the data. In addition, Cyg A and Cas A were modelled 
and subtracted per frequency channel to properly deal with the 
frequency dependence of the instrumental polarization. 

Since most of the Galactic polarized emission appears on 
spatial scales of a few arcmin, we make low resolution polar- 
ized images by discarding the antennas 0, C and D. With this 
selection the aiTay has ~4.2 arcmin resolution, where the con- 
version factor at this resolution is 1 mJy beam"' = 1 K. 

For the output RM cube, we selected the RM interval 
[-100, 100] rad m"^ where the Galactic emission most likely 
appears. From the RM synthesis cube we found that all the po- 
larized emission is concentrated in the interval [-7,1] rad m'^, 
which means that the line-of-sight component of the Galactic 
magnetic field is mainly pointed away from the observer. 
Figure [8] shows the most significant frames of the RM cube. 
Most of the frames still show unsubtracted polarized point 
sources at the level of a few mJy. In particular they appear to 
be present at the edge of the field, where the instrumental po- 
larization is stronger. 

The polarized emission shows a very peculiar trend as a 
function of Faraday depth. At -5 rad m"^ the emission is lo- 
cated in the center of the image, describing a bubble with a di- 
ameter of ~ 2°. The frame at -3 rad m"^ looks completely dif- 
ferent: in the centre of the image, at the location of the peak of 
the intensity, there is now a hole in the emission, suiTounded al- 
most completely by diffuse polarization. This pattern becomes 
even more evident at -2 rad m"^ before disappearing again at 
RM = 0. This frame shows the peak of the polarized emission 
together with the greatest spread throughout the field. There 
is an evident orientation in the spatial pattern with 'stripes' of 
emission running from South- West to North-East almost par- 
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allel to lines of constant Galactic latitude. This pattern was al- 
ready apparent in the 350 MHz data in HKB. 

In Figure|9]we show a few examples of polarized emission 
as a function of Faraday depth which further show the complex- 
ity of the interstellar medium through various lines of sight. 
In the first panel a bright point source lies at RM - 0. The 
slight asymmetry of the first side lobe might already indicate 
that the source is slightly more complicated than it looks. The 
second panel shows a double source: two peaks of the same 
intensity appear at -5 rad m"^ and 1 rad nT^. The third panel 
shows an even more complicated profile where multiple peaks 
are present at different RM values. 

It is possible to compute the total polarized intensity P by 
integrating along the Faraday depth. According to Brentjens 
(2007), it can be computed as follows: 

where crp ~ 1.4 K is the noise in polarization, P, is the po- 
larized map at the RM value / and the sum is over the frames 
of the RM cube which contain relevant emission. The factor B 
represents the area of the restoring beam divided by the inter- 
val between two frames of the RM cube. In our case we do not 
deconvolve the RM cube so we use the Gaussian which best 
fits the peak of the dirty beam ignoring the side lobe contribu- 
tion. This approximation is justified by the fact that most of the 
polarized flux in the individual frames of the RM cube is about 
7 K or lower, so its side lobe noise is lower than the thermal 
noise. Figure[TO]shows the integrated polarized intensity. 

Considering the inner 6° x6° square where we corrected for 
the primary beam shape, the nns of polarization fluctuations is 
7.2 + 0.8 K. 

4.1. Comparison with tlie 350 MHz WSRT results and 
the nature of the 'ring' 

We compared our polarization results with the 350 MHz data 
which have similar angular resolution. The most striking fea- 
ture present in the 350 MHz data is a ring-like structure which 
was first detected at 408 MHz in the late sixties (Bingham & 
Shakeshaft 1967) and then observed again at 21 cm with the 
Green Bank telescope (Verschuur 1969). 

In the 350 MHz data, the ring is mostly visible in polariza- 
tion angle. HKB did not perform RM synthesis but determined 
the rotation measures via the traditional linear fit to A^. They 
found that the RM values corresponding to the ring are mostly 
in the interval -8 < RM < -4 rad m"^ and find positive RM 
values - up to 12 rad nT^ - in the North-East and South-West 
corners of the field. 

In their work they concluded that the ring is probably a 
magnetic feature in the field parallel to the line of sight. They 
found no evidence of total intensity emission associated with 
the detected polarized structure. 

This ring-like structure is still present at 2 m wavelength, 
centred at o- ~ 3*' 15™, 6 ~ 65° and seen in emission for 
RM < -4 rad m"^ and as a lack of emission for -3 < RM 
< -1 rad m"^. 
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Fig. 9. Examples of Faraday spectra along different lines of 
sight together with the RMSF (bottom figure) of our observa- 
tions. 



12 



G. Bernardi et al.: WSRT observation of Galactic emission at 150 MHz 



70' 



68' 



o 






o 






o 


66° 




02 




H 


d 






o 






• i-H 












d 
•1 — 1 


64° 




"o 






0) 






Q 








62° 





60' 



7^ — \ — r 




J L 




- - 



4 



cd 
a; 

m 



1-^ 



3^30™ 



3^15™ 



3hQQm 



2'^45"^ 



Right Ascension (J2000) 

Fig. 10. Image of the total polarized intensity integrated along the Faraday depth. The conversion factor is 1 mJy beam ' = 1 K. 



The map of polarization angles shown in Figure [TT] makes 
the ring more evident. We can observe two regimes in the spa- 
tial distribution of polarization angles. Outside the ring, polar- 
ization angles are uniform on scales of several tens of arcmin, 
whereas there are variations in the polarization angle on much 
smaller scales inside it. At 350 MHz instead, polarization an- 
gles show structures on bigger scales inside the ring. 

The structure that we observe as a ring could therefore be 
due to a change in the spatial gradient in the polarization an- 
gle and this gradient could also be responsible for depolarizing 
the signal inside the ring at -3 < RM < -1 rad vnT^ where the 
polarized emission disappears. For the rest, the large scale dis- 



tribution in polarization angles displays the same orientation as 
the polarized intensity, from South-West to North-East. 

There are two differences between the 150 MHz and the 
350 MHz data. The 2 m data show almost no positive RM val- 
ues. This discrepancy might be due to the fact that the 350 MHz 
data are mosaic observations and, therefore, could lack the 
large scale emission, although HKB reject this possibility. A 
missing large scale structure can generate changes in the gradi- 
ent of the spatial distribution of RMs (Schnitzeler et al. 2007). 

The second important difference is the total intensity coun- 
terpart. Although there is no direct and strong correlation be- 
tween total intensity and polarization, in Section [3] we showed 
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that there are fluctuations in the Galactic radio background on 
scales between 13 arcmin and one degree. Assuming a spectral 
index (3-2.5 for the synchrotron emission, these fluctuations 
scale down to ~ 1.7 K at 350 MHz. This value is still a factor 
~2 above their noise and should allow a partial detection. The 
cause of this discrepancy between the 150 MHz and 350 MHz 
data remains unclear 

In light of the picture we have outlined so far, we suggest 
an origin for the nature of the ring different from the explana- 
tion proposed by HKB. They interpreted the ring as a magnetic 
structure, like a tube, which generates the variations of the mag- 
netic field along the line of sight responsible for the variations 
in the observed RM. 

At 150 MHz we observe diffuse total intensity emission 
and, therefore, we detect fluctuations in the perpendicular com- 
ponent of the magnetic field, responsible for the synchrotron 
emissivity. We do not, however, observe a tight correlation be- 
tween the polarized and total intensity components. It is hard 
to imagine a mechanism which generates fluctuations only in 
the parallel component of the magnetic field, leaving the per- 
pendicular one unperturbed. 

What is easier to imagine is a picture where at least two 
screens are acting in front of background synchrotron emis- 
sion. The background emission is Faraday rotated by a large 
scale screen which originates the bulk of the emission at RM = 
[-2, 1] rad m"^. In front of this, another perturbation in the dis- 
tribution of the thermal electrons along the line of sight, like a 
'bubble', would further rotate the polarized intensity. FigurefT2l 
sketches this scenario. 

This possible explanation would be supported by the mor- 
phology observed in the frames of the RM cube. As we men- 
tioned above, the frames at -3 rad m"^ and -2 rad m"^ show 
a cavity where the frame at -5 rad m"^ shows emission. If a 
foreground bubble with higher electron density was present in 
front of polarized background emission, it could just shift the 
polarized emission to higher RM values, as we observe. 

This only represents a qualitative interpretation of the data, 
and quantitative simulations that can predict the observed dis- 
tribution of polarized emission as a function of Faraday depth 
in the interstellar medium are needed in order to obtain a more 
constraining picture. 

5. Power spectrum analysis 

In Section [3] and Section |4] we have described fluctuations in 
the diffuse emission through the rms at a certain angular scale. 
We also argued that the total intensity fluctuations at arcmin 
scales are affected by source confusion and, therefore, not fully 
representative of the diffuse emission alone. In this section we 
present a power spectrum analysis in order to characterize com- 
pletely the emission as a function of the angular scale. 

The angular power spectrum is commonly used in cosmol- 
ogy to define the spatial properties of diffuse radiation. In in- 
terferometric observations, the power spectrum is usually com- 
puted in the uv plane, because the interferometer directly mea- 
sures the Fourier components of the brightness temperature of 
the sky (White et al. 1999; Myers et al. 2003; Morales & Hewitt 
2004). 




background emission 
of intrin.sic synclirotron 
polarized radiation 

Fig. 12. A possible picture of the line-of-sight magnetized in- 
terstellar medium in the Fan area. The intrinsically polarized 
background emission passes through multiple Faraday screens 
which make the spatial distribution of the polarized emission 
more clumpy. These screens generate the bulk of the polarized 
emission at RM = [-2, 1] rad mT^. A bubble of ionized ther- 
mal plasma sits between the Faraday screens and the observer. 
This bubble further rotates the polarized emission pushing it to 
RM —5 rad m l 

However, in this paper we compute the power spectra from 
the image as this approach allows for two simplifications. First, 
the correction for the primary beam of the interferometer can 
be applied much more easily in the image plane than in the uv 
plane. In the image plane the primary beam represents a sim- 
ple multiplicative effect whereas it is a convolution in the uv 
plane (Thompson et al. 1986). Second, the total polarization 
integrated along the whole Faraday depth is more easily com- 
puted in the image plane. 

Given this, we estimated the angular power spectrum as 
(Seljak 1997): 

" I ^ Z ^^'^^'^'^ - ^"'^^^ 

where X indicates either the total intensity / or the polarized 
intensity P,{ - ^ where is the angular scale in degrees, f2 
is the solid angle in radians, Nc is the number of Fourier modes 
around a certain ( value, X and X* are the Fourier transform of 
the image and its complex conjugate respectively, I is the two 
dimensional coordinate in Fourier space, cr^o^^^ is the rms noise 
and Nb is the number of independent synthesized beams in the 
map. The factor b^{€) is the power spectrum of the window 
function (Tegmark 1997). Since interferometric images repre- 
sent the true sky brightness convolved with the dirty beam (the 
Fourier transform of the weighted uv coverage), in our case the 
factor b^{C) is the power spectrum of the dirty beam. 

The number of modes around a certain ( value depends on 
the bin width in Fourier space and has a minimum dictated by 
the width of the field of view of the instrument. The 25 m dish 
of the WSRT telescope gives ~ 40, therefore we choose 
M - 50 as the bin width for computing the power spectrum. 
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Fig. 11. Images of polarization angles for the frames at RM = -2 rad nT^ (left) and RM = rad m"^ (right). In the color scale, 
black corresponds to -90° and white to +90°. 



The rms noise can be measured directly from the map, as 
we did in Section[3] However, since we have corrected the map 
for the primary beam, the noise was affected by the correc- 
tion as well. In order to understand how the expected noise 
power spectrum changed under the primary beam correction, 
we generated 100 realizations of a 6° x 6° noise map where 
the pixel values were drawn from a Gaussian distribution with 
o"noise - 0.75 mjy beam"' on 2 arcmin scales. Afterwards we 
multiplied the map by the primary beam shape. We computed 
a power spectrum for each noise map, and for each noise map 
after correcting for the primary beam shape. The resulting av- 
erage power spectra are shown in Figure [T3] 

We can see that the power spectrum remains flat at all the 
angular scales even after the map is multiplied by the primary 
beam function, but the amplitude increases. We find that the 
ratio between the two spectra is ~3.7, therefore we account for 
this factor in subtracting the rms noise in Equation|5] 

Finally, it is relevant to note that Equation |5] relates the 
power spectrum and the rms fluctuations T,^,, measured from 
a map as: 



In 



(6) 



where A^bin is the number of bins used to compute the power 
spectrum. 

Figure [14] presents the power spectrum down to 3 arcmin, 
where we avoided the smallest scales which are the most af- 
fected by residual calibration and ionospheric errors as we 
noted in Section[3] 

The total intensity power spectrum shows two clearly dif- 
ferent behaviours as a function of the angular scale. At large 
angular scales there is a power-law behaviour, with power de- 



Fig. 13. Solid line: average power spectrum computed on a 
6° X 6° noise map (see text for details). Triple-dot-dashed line: 
expected power spectrum of the thermal noise. Dashed line: 
average power spectrum of the noise after the primary beam 
correction is applied. 

creasing at increasing angular scale. This behaviour is typical 
of the Galactic diffuse emission observed at higher frequencies 
and higher angular resolutions (Bennett et al. 2003, La Porta et 
al. 2008). 

Around i ~ 900 the power spectrum flattens and remains 
flat down to 3 arcmin scales. This behaviour is expected if point 
sources dominate the emission. 

We have found in Section|3]that the image made by discard- 
ing the short baselines contains no diffuse emission, therefore 
we expect the power spectrum to be dominated by point sources 
at high t values. 
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Fig. 14. Filled circles: total intensity power spectrum from the 
residual image of Figure |2] with the best power-law fit super- 
imposed (solid line). Asterisks: polarization power spectrum 
from the total polarized intensity map of Figure [TO] with the 
best power-law fit superimposed (solid line). Dot-dashed line: 
power spectrum due to residual point sources (see text for de- 
tails). The plotted 1 cr error bars only account for the statistical 
errors. Power spectra are computed in the inner 6° x 6° square 
of the map. 

The angular power spectrum due to unsubtracted point 
sources is expected to be flat if the point sources have a ran- 
dom spatial distribution (Tegmark & Efstathiou 1996), there- 
fore it can be modelled as an additional noise term once its rms 
is known. We used the cr^^ value measured from the map of 
FigurelHas representative of the rms fluctuations due to unsub- 
tracted point sources and estimated the corresponding angular 
power spectrum as: 



(7) 



In Figure [14] the power spectrum expected from residual 
point sources agrees very well with the power spectrum com- 
puted from the map for { > 900, corresponding to ~12 arcmin. 
We can therefore conclude that the power spectrum down to 
{ ~ 900 is representative of the Galactic diffuse foreground 
emission whereas it is dominated by residual point sources at 
higher £ values. The sensitivity to diffuse structure on small an- 
gular scales is therefore not limited by instrumental noise, but 
by source confusion. The analysis presented in Section [3] fully 
agrees with the power spectrum analysis. 

The power spectrum of the total intensity diffuse emission 
was fitted by a power-law down to ^ = 900: 



•^f- '^400 ^55 



(8) 



The best fit values are C'^^ 



= 0.0019 + 0.0003 K2 and = 
-2.2 ± 0.3. The slope is within the range of slopes [-2,-3] 
found by all the previous measurements at higher frequencies 
(Tegmark & Efstathiou 1996, Giardino et al. 2001, Bennett et 
al. 2003, La Porta et al. 2008). 



The polarization power spectrum is fainter than the total in- 
tensity one and shows a regular decrease from the largest scales 
down to 4 arcmin scales. The power spectrum analysis con- 
firms that polarization data are not limited by confusion noise 
because of the very few polarized point sources, as the RM syn- 
thesis analysis already showed. The rms value down to 4 ar- 
cmin therefore represents fluctuations in the diffuse Galactic 
polarized foreground. 

Also the polarization power spectrum was fitted by a 
power-law down to ^ = 2700: 



{ 
700 



(9) 



The best fit values are C^o„ = 90 + 7 (mK)^ and jSf = 
-1.65 + 0.15. We note that the slope of the polarization power 
spectrum is flatter than the slope of the total intensity as was 
found by most of the studies of the synchrotron emission as a 
foreground of the cosmic microwave background (Tucci et al. 
2000; Baccigalupi et al 2001; Bruscoli et al. 2002; Tucci et al. 
2002). 

6. Discussion and conclusions 

We have presented results from Westerbork observations of the 
Fan region at 150 MHz both in total intensity and in polariza- 
tion, mainly focused on studying the foregrounds for the cos- 
mological 21 cm line. Our observations are the deepest avail- 
able so far, reaching a thermal noise of 0.75 mJy beam"'. 
However, on arcmin scales the main limitation comes from the 
classical confusion noise, which is ~3 mJy beam"' in the map. 

For the first time total intensity and polarization fluctua- 
tions in the Galactic diffuse foreground emission have been 
measured at the frequencies and angular scales relevant for 
EoR experiments. 

We have detected structure in the diffuse total intensity 
Galactic foreground at 150 MHz with an rms fluctuation of 
14 K at 13 arcmin resolution. The power spectrum analysis 
showed that the signal follows a power-law behaviour on the 
largest scales down io ( ~ 900 where it becomes compatible 
with the expected contribution from residual point sources. 

The best power-law fit to the power spectrum of Galactic 
emission gave an amplitude of C40Q - 0.0019 + 0.0003 K? at 
( = 400 and a slope I3\ = -2.2 ± 0.3. 

We have used the RM synthesis technique to measure the 
polarized emission and we found that polarization is quite com- 
plex and structured both spatially and in Faraday depth. The 
RM synthesis technique has confirmed its status as a superb 
tool to detect multiple layers of emission at different Faraday 
depths. 

Although the large scale distribution of polarized emission 
somewhat resembles the total intensity one, the lack of corre- 
lation on intermediate and small angular scales suggests that 
the polarized emission is generated through the Faraday screen 
mechanism. 

The 'ring' already observed at 350 and 408 MHz has been 
detected again in our data. From the information available at 
different Faraday depths, we suggest that the ring is more likely 
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to be a feature in the density of thermal electrons, like a fore- 
ground bubble, rather than a magnetic structure. 

The detection of polarization at 150 MHz indicates that the 
local medium must be Faraday thin, otherwise most of the po- 
larization structures along the line of sight would vanish at 2 m 
due to Faraday depolarization. Some Faraday depolarization is 
also present in the Fan area, however, because we would expect 
~40 K of polarized signal on the basis of a direct extrapolation 
of the 350 MHz data. Instead we found 7.2 K of nns fluctua- 
tions at 4 arcmin after integrating along the Faraday depth. 

The polarization angular power spectrum shows a mono- 
tonic decrease from the largest angular scales down to 4 ar- 
cmin. All the polarized signal comes from diffuse structure in 
the Galaxy. The best power-law fit to the power spectrum of 
polarized emission gave an amplitude of C^qq = 90 + 7 (mK)^ 
at ^' = 700 and a slope ySf = -1.65 ± 0.15. 

The fact that the Fan region shows a strongly polarized sig- 
nal over scales ranging from several degrees down to 4 arcmin 
scales makes it a promising target to study the behaviour of the 
ionospheric Faraday rotation and to calibrate the instrumental 
polarization of the forthcoming low frequency arrays. 

The results presented so far can be used to improve the pre- 
dictions of the level of foreground contamination of the EoR 
signal. 

The Fan region is obviously not a potential target of EoR 
observations, because of its proximity to the Galactic plane. 
However, the power spectrum of total intensity fluctuations 
measured there can be extrapolated down to smaller angular 
scales and used to obtain an indication of the level of the fluc- 
tuations at higher Galactic latitude, where EoR observations 
will certainly be carried out. We will limit ourselves to present- 
ing results in terms of two-dimensional power spectra which 
only describe the spatial properties of the sky brightness. The 
EoR signal has an intrinsic three dimensional nature which will 
ultimately be taken into account in order to remove the fore- 
grounds and to characterize the statistics of the cosmological 
signal (Zaldarriaga et al. 2004, Datta et al. 2007, lliev et al. 
2008), but we postpone this analysis to future work. 

In all the radio surveys of the Galaxy, the diffuse emis- 
sion decreases away from the Galactic plane (see de Oliveira- 
Costa et al. 2008 for a recent analysis of multi-frequency all- 
sky data). In Figure|6]we plotted the behaviour of the nns fluc- 
tuations as a function of Galactic latitude in the 2° resolution 
all-sky map at 150 MHz and showed that the signal drops by 
a factor of ten already at moderate Galactic latitude, staying 
roughly constant up to the Galactic pole. Also in surveys at 
higher frequencies and with slightly better angular resolution 
(Jonas, Baart & Nicolson 1998; Gold et al. 2008), it is possible 
to find "cold" spots where the signal is ten times fainter than in 
the proximity of the plane or even lower 

Therefore, if we extrapolate the measured power spectrum 
of diffuse total intensity emission through its best power-law 
fit down to 5 arcmin - around where the EoR signal is ex- 
pected to peak - we find that its level at ^ = 2160 is (jr = 

+ \)C'^l2n ~ 5.7 K and the corresponding rms value is 
~ 1 8.3 K. At moderate and high Galactic latitudes this value can 
drop by a factor of ten or even more if favourable sky patches 
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Fig. 15. Comparison between angular power spectra of various 
components relevant for the EoR detection. From top to bot- 
tom: total intensity (solid line) and polarization (short-dashed 
line) power spectra fitted from the Fan data and extrapolated 
down to ~3 arcmin, power spectrum due to unresolved point 
sources (long-dashed line), extrapolation at high Galactic lati- 
tudes of the best fit power spectrum of the total intensity emis- 
sion (dot-dashed line), power spectrum of the thermal noise 
of our observations (triple-dot-dashed line), theoretical expec- 
tations for the EoR signal (Zaldarriaga et al. 2004, downward 
arrows). 



are identified. Figure [T5]compares the measured power spectra 
of foregrounds and the expected EoR signal. 

It is worth noting in this context that the detection of diffuse 
foregrounds is already limited by source confusion at arcmin 
scales in an area like the Fan region, where the diffuse emis- 
sion is quite bright. This indicates that source confusion could 
be the most serious foreground contaminant in detecting the 
cosmological signal at arcmin scales. In Figure [TSl we plotted 
the power spectrum due to residual point sources assuming that 
its normalization is given by the rms value of ~3 mJy beam"' 
found in Section [3] This contribution might dominate com- 
pletely over the diffuse emission at angular scales smaller than 
~5 arcmin at all Galactic latitudes. 

Future experiments aimed at detecting the EoR signal 
would therefore benefit from the inclusion of long baselines 
in order to decrease the confusion limit. 

Finally, we have already argued that the polarization sig- 
nal has also to be considered in the context of the EoR detec- 
tion, because a polarized signal coming from the sky will leak 
into Stokes /, inducing artificial fluctuations on top of the cos- 
mological ones if the calibration is not perfect. Moreover, the 
leaked signal would vary with frequency because of Faraday 
rotation, generating a signal which would have characteristics 
similar to the cosmological one and, therefore, be extremely 
difficult to remove. 

The seriousness of this phenomenon depends upon the 
strength of the polarized signal coming from the sky, its RM 
- higher RM values will induce more rapid variations across 
the bandwidth - and the calibration accuracy. 
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Recent observations at 150 MHz with the GMRT have 
found a region of the sky 1 ° wide where the polarized bright- 
ness temperature is less than 1 K down to ^ ~ 1000 (Pen et al. 
2008), suggesting that sky regions free of synchrotron polariza- 
tion exist. We note however, that they do not impose constraints 
on 5-10 arcmin scales. 

We have detected nns fluctuations in the difiiise polarized 
emission of 7.2 K at 4 arcmin. On 5 arcmin scales the level of 
the polarization power spectrum is ST ~ 3.3 K. Most of the 
polarized emission appears at RM values [-7,1] rad m"^. Such 
values do not generate rapidly rotating Stokes Q and U signals 
over a ~2 MHz bandwidth but they can potentially mimic an 
EoR signal over a 5-7 MHz bandwidth. 

Generalizing this result appears to be rather difficult. The 
Fan region is certainly a peculiar region in polarization. Its po- 
larization power spectrum cannot be taken as representative of 
other sky areas since the region centred around the ring con- 
tributes most to it. 

It is reasonable to beUeve that the polarization power spec- 
trum represents an upper limit on most of the sky at moderate 
and high Galactic latitudes, where lower values should be ex- 
pected, in agreement with what the preliminary GMRT results 
suggested. 

It is worth pointing out, however, that whereas the polar- 
ized brightness temperature in the Galactic halo should be in- 
trinsically lower than in the Fan area, the interstellar medium 
is probably Faraday thin, because of lower values of the elec- 
tron density and magnetic field. Therefore regions outside the 
plane could also be structured in Faraday depth and RM syn- 
thesis will be the way to detect potentially weak and structured 
polarized signals. 

Future work will be dedicated to further characterizing 
foregrounds at higher Galactic latitude, with particular atten- 
tion paid to their polarization features. 
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